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Quantum Monte Carlo calculations of the relaxation energy, pair-correlation function, and 
annihilating-pair momentum density are presented for a positron immersed in a homogeneous elec- 
tron gas. We find smaller relaxation energies and contact pair-correlation functions in the important 
low-density regime than predicted by earlier studies. Our annihilating-pair momentum densities 
have almost zero weight above the Fermi momentum due to the cancellation of electron-electron 
and electron-positron correlation efi'ects. 
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Electron-positron annihilation underlies both medical 
imaging with positron emission tomography (PET) and 
studies of materials using positron annihilation spec- 
troscopy (PAS) [ll. Positrons entering a material rapidly 
thermalize and the majority annihilate with opposite- 
spin electrons to yield pairs of photons at energies close 
to 0.511 MeV. In a PET scan, positrons arc emitted by 
radionuclides in biologically active tracer molecules and 
the resulting annihilation radiation is measured to image 
the tracer concentration. The interaction of low-energy 
positrons with molecules is therefore of substantial exper- 
imental and theoretical interest [Ij. PAS is used to in- 
vestigate microstructures in metals, alloys, semiconduc- 
tors, insulators [l|, polymers and nanoporous mate- 
rials Q . Positrons are repelled by the positively charged 
nuclei and tend to become trapped in voids within the 
material. The positron lifetime is measured as the inter- 
val between the detection of a photon emitted in the /3+ 
radioactive decay that produces the positron and the de- 
tection of the annihilation radiation [ij . The lifetime is 
characteristic of the region in which the positron settles, 
and PAS is a sensitive, nondestructive technique for char- 
acterizing the size, location, and concentration of voids 
in materials. Measuring the Doppler broadening of the 
annihilation radiation or the angular correlation between 
the two 0.511 MeV photons yields information about the 
momentum density (MD) of the electrons in the pres- 
ence of the positron. These techniques may be used to 
investigate the Fermi surfaces of metals [Hj . 

The aim of PAS experiments is to investigate a host 
material without the changes induced by the positron. 
The positron is, however, an invasive probe which polar- 
izes the electronic states of the material. Disentangling 
the properties of the host from the changes induced by 



the positron is a major theoretical challenge. Positrons in 
condensed matter may be modeled with two-component 
density functional theory (DFT) Q, in which the cor- 
relations arc described by a functional of the electron 
and positron density components. Within the local den- 
sity approximation (LDA), this functional is obtained 
from the difference Afi between the energy of a homoge- 
neous electron gas (HEG) with and without an immersed 
positron. Afi is known as the relaxation energy, and is 
equal to the electron-positron correlation energy. 

Two-component DFT gives reasonable electron and 
positron densities, but the DFT orbitals do not de- 
scribe electron-positron correlation properly [1,0. The 
electron-positron pair-correlation function (PCF) g{r) 
and the annihilating-pair momentum density (APMD) 
p{p) constructed from the DFT orbitals are therefore 
poor. The contact PCF 17(0) is particularly impor- 
tant because it determines the annihilation rate A = 
35(0)/(4c^?'s) [H for a positron immersed in a param- 
agnetic HEG, where Ts is the electron density parameter 
and c is the speed of light in vacuo 0] . If the electron and 
positron motions were uncorrclatcd g{Q) would be unity, 
but the strong correlation leads to much larger values, 
particularly at low densities, where an electron-positron 
bound state (positronium or Ps) or even an electron- 
electron-positron bound state (Ps~) may be formed. 

We have used the variational and diffusion quantum 
Monte Carlo (VMC and DMC) methods 0, [ifl as im- 
plemented in the CASINO code [ll[ to study a single 
positron in a HEG. Fermionic antisymmetry is imposed 
via the fixed-node approximation, in which the nodal sur- 
face is constrained to equal that of a trial wave function. 
We used Slater- Jastrow (SJ) and Slater- Jastrow-backflow 
(SJB) trial wave functions [12, ll3|- The latter go be- 
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yond the single-particle SJ nodal surface by replacing 
the particle coordinates in the Slater determinants by 
"quasiparticle coordinates." SJB wave functions giv e the 
highest accuracy obtained to date for the HEG [13] ■ 
We also tested two types of orbitals: (i) plane-wave or- 
bitals for each particle and (ii) orbitals which describe the 
pairing between the electrons and positron. The pairing 
orbitals were obtained from mean-field calculations per- 
formed in the reference frame of the positron, so the or- 
bitals are functions of the separation of an electron and 
the positron [l4|. Within this impurity- frame DFT (IF- 
DFT) method, the pairing orbitals describe the electron- 
positron correlation quite well on their own IJ] and give 
a different nodal surface from the plane-wave orbitals. 
(NB, our QMC calculations were performed in the labo- 
ratory frame.) The four wave- function forms used are: 
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where R denotes the positions of all the particles, r-|- and 
denote the positions of up- and down-spin electrons, 
respectively, rp is the positron position, and [• • • ] de- 
notes a Slater determinant. The Jastrow exponent J(R) 
ISf and the backflow displacement ^(R) ll3| contain pa- 
rameters that were optimized separately for each wave 
function and system. The Jastrow exponents were first 
optimized usin g th e efficient VMC variance-minimization 
scheme of Ref. Il6|, and then all the parameters (includ- 
ing the backflow parameters) were optimized together 
using the VMC energy-minimization scheme of Ref. [TtI . 
The pairing orbitals {4>i} were represented using B-spline 
functions on a real-space grid [iSj . The electron-positron 
cusp condition was enforced on the pairing orbitals for 
wave function ^^l^^ [11, for the other three wave 
functions, the cusp conditions were imposed via the Jas- 
trow factor. In all our calculations the simulation-cell 
Bloch vector [2lj was chosen to be = 0. 

Tests at high (r^ = 1) and low (r^ = 8) electron den- 
sities show that the qualitative features of the variations 
in AQ, g(r), and p(p) with Ts are the same for each of 
the four wave functions of Eq. ([1]). However, as shown in 
the auxiliary material [22|, we obtained lower VMC and 
DMC energies with the SJB wave functions ('I'l'^ ^^'^^ 



^'^ ), and therefore 



^f^^g) than the SJ ones (^'pw and ^'p^^^^ 
we used SJB wave functions to obtain all our main re- 
sults. The pairing orbitals give lower SJB- VMC energies 
than the plane-wave orbitals, but the SJB-DMC ener- 
gies with the plane-wave and pairing orbitals arc almost 
identical. The lack of sensitivity to the orbitals used, 
and hence the nodal surface, suggests that the DMC en- 
ergies are highly accurate. The energies reported in this 



paper are from DMC calculations using wave function 
^'l'^^. Such calculations are considerably less expensive 
than calculations using 'ip^ due to (i) the lower energy 
variance achieved with [l^] and (ii) the fact that 

plane-wave orbitals are cheaper to evaluate. The DMC 
energies were extrapolated to zero time step. Our pro- 
duction DMC calculations were performed in cells con- 
taining iV = 54 electrons. Tests of convergence with 
respect to system size up to iV = 114 electrons are de- 
scribed in the auxiliary material [2^. The cell volume 
was chosen to be {N — 1) x (4/3)7rrg, so that the elec- 
tron density far from the positron was correct. IF-DFT 
calculations 14| suggest that finite-size effects due to the 
interaction of images of the positron are negligible for 
> 54 electrons. 

Our DMC relaxation energies are plotted in Fig.[T]and 
are well-fitted by the form 



Ao + Airs - 0.262005B2r2 
1 + BiTs + B2r1 



(2) 



where A^i = -0.260361, Aq = -0.261762, Ai = 
0.00375534, Bi ^ 0.113718, and B2 = 0.0270912. Equa- 
tion ([2]) tends to the correct low-density limit of the en- 
ergy of the Ps~ ion . Equation ([2]) does not yield the 
exact high-density behavior of the random phase approx- 
imation (RPA), although this is only relevant for rs < 0.1 
[i^ l. VMC energies for a positron in a HEG have been 
reported previously [25| , but we have used superior trial 
wave functions and have obtained very different results. 
At high densities our relaxation energies are similar to 
those of Lantto [2^, but at lower densities we obtain 
smaller values. The SJB-DMC and IF-DFT results [ll 



and the data of Borohski and Stachowiak [27| show simi- 
lar behavior with Tg , while the Borohski- Nieminen fit [^ 
to the data of Ref. |28| is markedly different. Borohski 
and Nieminen's ^ expression for Ari(rs) is widely used 
in two-component DFT calculations, but our study sug- 
gests it is not very accurate and should be replaced by 
Eq. ©. 

We calculated the APMD within VMC using optimized 
SJB trial wave functions with pairing orbitals (v^pai^); 
because these give lower VMC energies than plane- wave 
orbitals (^'p'^^). These calculations were performed by 
constraining an electron and the positron to lie on top 
of one another throughout the simulation [2^. APMDs 
at different densities are plotted in Fig. [2l with the nor- 
malization chosen such that Jp 47rp^p(p) dp = (4/3)7rfc|n. 
Our results clearly show the enhancement of the APMD 
below the Fermi momentum predicted by Kahana [30j . 
but our data differ quantitatively from previous results 
14, 27, 11]. Our VMC data have almost no weight above 



the Fermi momentum over the entire density range stud- 
ied, even though the weight in the MD above /cf in the 
HEG is substantial at low densities. For example, we find 
that the APMD immediately above kp is roughly 10% of 
the value for the HEG at Tc = 1 and 3% at Tc = 8. 
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FIG. 1: (Color online) Relaxation energy against density pa- 
rameter from our SJB-DMC calculations and other studies 
[n.[26l.j28l. [29l|. relative to the Bororiski-Nieminen expression, 
AQbn [g| (horizontal dashed line). 



Suppression of the weight in the APJMD above kp was 
demonstrated theoretically by Carbotte and Kahana [sij , 
but our study gives a more detailed and accurate pic- 
ture. We investigated the weight above kp using VMC 
with the wave function ^'p'^^ by selectively eliminating in- 
terparticle correlations. Neglecting electron-electron and 
electron-positron correlations gives the familiar "top hat" 
JMD of the noninteracting system. Calculations with the 
electron-positron terms removed give an APMD indis- 
tinguishable from the MD of the HEG, with a tail above 
kp. Calculations including electron-positron correlation 
but neglecting electron-electron correlation show Kahana 
enhancement below kp and a tail above kp. When, how- 
ever, both electron-electron and electron-positron corre- 
lations are included, the tail above kp is largely sup- 
pressed, as shown in the lower panel of Fig. [2l 

The suppression of the tail in the APJVID can be ex- 
plained by examining the behavior of the two-body terms 
in the Jastrow exponent. (For simplicity, we consider the 
^pw wave function in the following discussion.) The 
Jastrow exponent J(R) is the sum of electron-electron 
[■ui-|-(r) and u^i{r), where the arrows indicate spins] and 
electron-positron [uop(r)] terms. If one assumes that 



(3) 



then the APMD has exactly zero weight above kp, as 



shown in the auxiliary material [22]. The RPA (linear 
response theory) shows that Eq. ^ holds at large r and 
the Kato cusp conditions force the gradients of u^iir) and 
Ucp{r) to satisfy Eq. ([3]) at r = 0. The cusp conditions for 
parallel and antiparallel spin electrons are different and 
therefore u^i{r) and M-t-t('^) must differ at small r, but an- 
tisymmetry ensures that the probability of parallel-spin 
electrons being closer than is small. As shown in the 
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FIG. 2: (Color online) Top: APlVIDs [p{p)] for different densi- 
ties. The solid lines show our VMC data obtained with wave 
function '^pi^ and = 114 electrons, while the dashed and 
dashed-dotted lines show the data of Kahana [sO] and Sta- 
chowiak [231, respectively. Bottom: APMDs for the positron- 
in-HEG and the HEG at = 8 and A'^ = 54 electrons, calcu- 
lated using *|w- 



auxiliary material, plots of the terms in the Jastrow ex- 
ponent demonstrate the approximate validity of Eq. . 

We calculated the PCFs within VMC and DMC using 
^pw wave functions, because these give the same results 
as pairing orbitals but the calculations are much cheaper 
[22j . The final results were evaluated by extrapolated 
estimation (twice the DMC PCF minus the VMC PCF) 
[3^ . in order to eliminate the leading-order errors. In 
Fig. ini the electron-positron contact PCF g(0) is plotted 
relative to the Boronski- Nieminen forrn [29i. which is a 
fit to the data of Stachowiak and Lach [33[ . Our contact 
PCF data are well-represented by 



.g(0) = l + 1.23r, ^ 

r.8/3 



+ aar^ + a^/gr 



■ 08/3? 



a3/2f, 

-0.173694r3 
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(4) 



where 



03/2 



-3.38208, a2 



.6957, 



0-7/3 



-7.37037, 
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and ttg/s = 1.75648. Equation ^ satisfies the high- 
density (RPA) [gj] and low-density (Ps~) hmiting be- 
haviors [23[. Our fuh data for g{r) are given in the aux- 
ihary material . The IF-DFT data follow the extrap- 
olated SJB data quite well, while the other many-body 
calculations give somewhat larger values of g{r) at low 
densities. In the density range = 5-8 a.u., our values 
of g{0) are approximately 9% smaller than those given 
by the Boronski-Nieminen expression The local in- 
crease of the electron density around the positron caused 
by their mutual attraction is modeled in two-component 
DFT using an "enhancement factor" based on data for 
g{0). Using our smaller values of ^(0) would reduce the 
enhancement factor and hence the overestimation of an- 
nihilation rates obtained with the positronic LDA jsj]. 
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FIG. 3: (Color online) Deviation of the contact PCF g{0) 
from the form pbn(O) of Boronski and Nieminen [29|| toRether 
with other results in the literature [3, [H, H^, l33l. |35|. [3q | . 



In conclusion, our results are the most accurate ob- 
tained so far for a positron in a HEG. Our data for 
AJ7 are sufficient to define the energy functional for a 
two-component positronic DFT within the LDA. They 
would also be useful in developing semilocal [s^l or other 
functionals. Our PCF data give a smaller enhance- 
ment factor than the standard Boronski-Nieminen ex- 
pression 0. Our APMDs have very little weight above 
kp because of the cancellation of electron-electron and 
electron-positron correlation effects. We have derived an 
exact result relating Eq. (jSj to the complete absence of 
weight in the APMD p(p) for p > kp, which is useful in 
understanding this effect. 

We acknowledge financial support from the UK En- 
gineering and Physical Sciences Research Council (EP- 
SRC). Computer resources were provided by the Cam- 
bridge High-Performance Computing Facility and the 
Lancaster High-End Computing cluster. 



[6: 

[r 

[8 
[9 

[lo: 
[11 

[12 

[13 

[14 

[15 

[16 

[17 

[18 
[19 

[20 

[21 
[22; 



[23' 
[24 
[25^ 



[26 
[27 
[28' 
[29; 

[30 
[31 

[32] 



R. Krause-Rehberg and H.S. Leipner, Positron Annihila- 
tion in Semiconductors (Springer- Verlag, Berlin, 1999). 
G.F. Gribakin, J. A. Yomig, and C.M. Smko, Rev. Mod. 
Phys. 82, 2557 (2010). 

R.A. Pethrick, Prog. Polym. Sci. 22, 1 (1997). 

D. W. Gidley, H.-G. Peng, and R.S. Vallery, Annu. Rev. 
Mater. Res. 36, 49 (2006). 

Zs. Major, S.B. Dugdale, R.J. Watts, G. Santi, M.A. 
Alam, S.M. Hayden, J.A. Duffy, J.W. Taylor, T. Jarl- 
borg, E. Bruno, D. Benea, and H. Ebert, Phys. Rev. 
Lett. 92, 107003 (2004). 

E. Boronski and R.M. Nieminen, Phys. Rev. B 34, 3820 
(1986). 

M.J. Puska and R.M. Nieminen, Rev. Mod. Phys. 66, 
841 (1994). 

We use Hartree atomic units {fi = |e| — — 47reo ~ 1) 
throughout. 

D.M. Ceperley and B.J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

W.M.C. Foulkes, L. Mitas, R.J. Needs, and G. Rajagopal, 
Rev. Mod. Phys. 73, 33 (2001). 

R.J. Needs, M.D. Towler, N.D. Drummond, and P. Lopez 
Ri'os, J. Phys.: Condens. Matter 22, 023201 (2010). 
Y. Kwon, D.M. Ceperley, and R.M. Martin, Phys. Rev. 
B 58, 6800 (1998). 

P. Lopez Ri'os, A. Ma, N.D. Drummond, M.D. Towler, 

and R.J. Needs, Phys. Rev. E 74, 066701 (2006). 

N.D. Drummond, P. Lopez Ri'os, C.J. Pickard, and R.J. 

Needs, Phys. Rev. B 82, 035107 (2010). 

N.D. Drummond, M.D. Towler, and R.J. Needs, Phys. 

Rev. B 70, 235119 (2004). 

N.D. Drummond and R.J. Needs Phys. Rev. B 72, 
085124 (2005). 

C. J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and 
R.G. Hennig, Phys. Rev. Lett. 98, 110201 (2007). 

D. Alfe and M.J. Gillan, Phys. Rev. B 70, 161101 (2004). 
A. Ma, M.D. Towler, N.D. Drummond, and R.J. Needs, 
J. Chem. Phys. 122, 224322 (2005). 

S.J. Binnie, S.J. Nolan, N.D. Drummond, D. Alfe, N.L. 
Allan, F.R. Manby, and M.J. Gillan, Phys. Rev. B 82, 
165431 (2010). 

G. Rajagopal, R.J. Needs, A. James, S.D. Kenny, and 

W.M.C. Foulkes, Phys. Rev. B 51, 10591 (1995). 

See EPAPS Document No. xxxxxxxxxx. 

For more information on EPAPS, see 

'http://www.aip.org/pubservs/epaps.htmlj 

A.M. Frolov, Phys. Lett. A 342, 430 (2005). 

J. Arponen, J. Phys. C 11, L739 (1978). 

G. Ortiz, PhD thesis, Swiss Federal Institute of Technol- 
ogy, Lausanne (1992). The relevant data from this thesis 
are also reported in Refs. [35l and[29l. 

L.J. Lantto, Phys. Rev. B 36, 5160 (1987). 

H. Stachowiak, Phys. Rev. B 41, 12522 (1990). 

J. Arponen and E. Pajanne, Ann. Phys. 121, 343 (1979). 

E. Boronski and H. Stachowiak, Phys. Rev. B 57, 6215 
(1998). 

S. Kahana, Phys. Rev. 129, 1622 (1963). 

J.P. Carbotte and S. Kahana Phys. Rev. 139, A213 

(1965). 

D.M. Ceperley and M.H. Kalos, in Monte Carlo Meth- 
ods in Statistical Physics 2nd edn, edited by K. Binder 



5 



(Springer- Verlag, Heidelberg, 1979), p. 145. 
[33] H. Stachowiak and J. Lach, Phys. Rev. B 48, 9828 (1993). 
[34] J. Mitroy and B. Barbiellini, Phys. Rev. B 65, 235103 

(2002). 

[35] V. Apaja, S. Denk, and E. Krotscheck, Phys. Rev. B 68, 
195118 (2003). 



[36] A. Harju, B. Barbiellini, S. Siljamaki, R.M. Nieminen, 
and G. Ortiz, J. Radioanal. Nucl. Chem. 211, 193 (1996). 

[37] B. Barbiellini, M.J. Puska, T. Korhonen, A. Harju, T. 
Torsti, and R.M. Nieminen, Phys. Rev. B 53, 16201 
(1996). 



